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Abstract 



We present a NLO study of inclusive polarized prompt photon production in a conceivable 
fixed target pp mode of HERA with longitudinally polarized protons at = 39 GeV. 
We analyze the sensitivity of the corresponding double spin asymmetry to the proton's 
polarized gluon distribution Ag and estimate the expected statistical precision in its 
determination. The main theoretical uncertainties in the predictions are examined. 



It has recently become possible to perform a complete and consistent study of longitu- 
dinally polarized deep-inelastic scattering (DIS) in next-to-leading order (NLO) of QCD, 
since the spin-dependent two-loop splitting functions, needed for the NLO evolution of 
the polarized parton distributions, have become available M, 0. A first such phenomeno- 
logical NLO analysis, taking into account all available experimental data on polarized 
DIS [[| has been presented in ||, followed by the analyses || ||. The studies of [§, |(| 
have shown that present polarized DIS data are still quite far from providing accurate 
knowledge about the nucleon's spin-dependent sea quark and gluon distributions. This 
holds true, in particular, for the polarized gluon density Ag, the x-shape of which seems 
to be hardly constrained at all [|], |6| by the DIS data, even though a tendency towards 
a sizeable positive total gluon polarization, L Ag(x, Q 2 = 4 GeV 2 )<ix > 1, was found 

[5], Thus, there is clearly some need for independent information on Ag. For this 
purpose, it seems expedient to look at processes for which Ag enters in leading order 
(LO) already, rather than as a NLO correction as for the spin-dependent DIS structure 
function g\. One of such processes is inclusive large-pr prompt photon production in 
collisions of longitudinally polarized protons, pp — > 7X j7|, |S], [CTJ. In the unpolarized 
case where this process has been studied in a huge number of experiments it has been 
an invaluable tool for pinning down the proton's unpolarized gluon distribution g(x, Q 2 ) 
]TT| , |T2| , |T3| , rj|. Hence prompt photon production with polarized beams seems a promising 



source for obtaining information on Ag. 

It is being discussed as one future option for HERA to polarize its 820 GeV proton 
beam. If this can be achieved, one could use the beam in a fixed target experiment, scat- 
tering it off an internal polarized nucleon target. This conceivable constellation, dubbed 
'Phase II' of HERA-N [T5|] , would yield y/s ~ 39 GeV and thus could provide information 



complementary to that obtained from planned similar spin physics experiments at much 
higher energies at the RHIC collider fl6|| . Theoretical predictions for polarized prompt 



photon production at y/s ~ 40 GeV have been made in the past || [K|, taking into ac- 
count the spin-dependent 'direct' subprocess cross sections for ab — > 7A (a, b = q,q,g) 
including their full NLO QCD corrections as calculated in || |9|. From the experience in 
the unpolarized case, the inclusion of NLO corrections is expected to be quite important 
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in order to make reliable predictions. The main shortcoming of the studies || [K| was, 
however, that spin-dependent parton densities evolved only in LO had to be used at that 
time. Having sets of NLO polarized parton distributions available now, we can obviously 
put the corresponding predictions as well as the assessment of their theoretical uncertain- 
ties such as their scale dependence on a much firmer basis, which is the purpose of this 
paper. The only remaining drawback here is that the fragmentation contribution to the 
polarized prompt photon cross section still cannot be calculated in NLO since the NLO 
corrections to the underlying polarized subprocesses ab — > cX (a, b,c = q, q, g) are not yet 
known. On the other hand, in the unpolarized case the fragmentation piece is known to be 



subdominant - even though not negligible - at fixed target energies (see, e.g., [14]). The 
fragmentation contribution to the polarized cross section which was omitted altogether in 
|| |TU| , will be included on a LO basis in this paper. Even though this is not strictly con- 



sistent in the framework of a NLO calculation, it is the best 'state-of-the-art' procedure 
and also appears reasonable in view of the fact that the K- factor K^ ra9 = &NLo/ a LO Q f° r 
the fully inclusive fragmentation part generally turns out to be very close to unity when 
calculated for the unpolarized case. 

The analyses of [|], []] provide several different sets of LO and NLO spin- dependent 
parton densities, all of which are in very good agreement with the existing polarized DIS 
data but differ mainly in the x-shapes of their polarized gluon distributions. We are 
therefore in the position to study the sensitivity of polarized prompt photon production 
to Ag(x,Q 2 ) on the basis of consistent NLO sets of parton distributions that include all 
experimental information presently available from DIS, but also reflect the full freedom 
concerning Ag left by those data. The latter is illustrated in Fig. 1 which compares 
the gluon distributions of various NLO sets of j|, || in the x-range dominantly probed 
by prompt photon production at a/s = 39 GeV and > 3 GeV. The sets we will 
use in our study are the NLO 'valence' set of the 'radiative parton model analysis' 0, 
which corresponds to the best-fit result of that paper (hereafter referred to as 'fitted Ag' 
scenario), and two other sets of 0] which are based on either assuming Ag(x, fi 2 ) = g(x, fi 2 ) 
or Ag(x,fi 2 ) = at the low input scale \i of ||, where g(x,fi 2 ) is the unpolarized NLO 
GRV [[HJ input gluon distribution. These two sets will be called 'Ag = g input' and 
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'Ag = input' scenarios, respectively, in what follows. It turns out that in the x-range 
explored here the NLO sets A,B of || have gluon distributions quite similar to those of 
the 'Ag = g input' and 'fitted Ag' scenarios of [|J, respectively. Only the gluon of set 
C of ('GS C) is qualitatively different since it has a substantial negative polarization 
at large x. We will therefore also use this set in our calculations. A graph similar to 
Fig. 1 could be shown for the corresponding LO gluon distributions of @, |6| (see, e.g., 
fT8fl). The fact that [|], |6| have provided both consistent LO and NLO parton sets in all 
cases, enables us to more reliably study the perturbative stability of the cross sections 
and asymmetries. This is again a major improvement with respect to the previous studies 
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The quantity to be studied in prompt photon experiments with polarized beam and 

target is the double-spin asymmetry 

^ _ da ++ — da + ~ _ dAa 
da ++ + da^ da ' 

where da ++ (da + ~) denotes the cross section for the prompt photon being produced by 

protons with same (opposite) helicities. Two types of processes contribute to the prompt 

photon production cross section: the so-called 'direct' piece, where the photon is emitted 

via a pointlike (direct) coupling to a quark, and the fragmentation piece, in which the 

photon originates from the fragmentation of a final state parton. The cross section for 

the fully inclusive production of a prompt photon with momentum p 1 thus schematically 

reads 

dAa = - {da ++ — da + ~) = dAa^ir + dAaf rag (2) 



^2 / dx a dx b Af a (x a , n 2 F )Af b (x b , n 2 F ) 

< 7Z 



a,b=q,q,g 



dAa2 b (p 7 , x a , x b , fj. R , li f , M f ) 



/dz 
— d ^° C ab(P~H X a, %b, Z, [Mr, fi F , M F )D2(z, M|) 

^=q,q,y 

where the dAa l ab represent the spin-dependent subprocess cross sections for partons a, h 
producing a particle i (i = 7, q, q, g), integrated over the full phase space of all other 
final state particles. The polarized parton distributions are defined as 

Afi(x, /4) = ff(x, fi 2 F ) - ff(x, ti F ) , (3) 



with f*(x, fj/p) (fi~(x, fip)) denoting the number density of parton-type i with momentum 
fraction x and positive (negative) helicity in a proton with positive helicity at scale y,p. 
Furthermore, in (@) D^(z, Mp) is the (unpolarized) photon fragmentation function at scale 
Mp, z being the fraction of energy of the fragmenting parton c transferred to the photon. 
We note that even in the polarized case the photon fragmentation functions are always the 
unpolarized ones since the polarization of the outgoing photon is not observed. Despite 
the fact that its corresponding partonic subprocesses are of order ct^, the fragmentation 
contribution is present already in LO since the parton-to-photon fragmentation functions 
are effectively of order a e , m Ja s in perturbative QCD. As made explicit in (|2|), the cross 
section in any fixed order of perturbation theory depends on unphysical scales which have 
to be introduced in the procedure of renormalization and of factorization of initial 
(Hf) and final (Mp) state mass singularities. Unless stated otherwise, we will choose 
Hr = Hf = Mp = p T /2 in what follows. 

The corresponding expressions for the unpolarized cross section, which we need to 
calculate the spin asymmetry A LL , can be obtained from by omitting the 'A' and 

taking sums instead of the differences on the rhs. In this case, all cross sections da l ab (i = 7, 
q, q, g) are known to NLO accuracy, the corrections to the direct and fragmentation 
subprocess cross sections having been calculated in the MS scheme in [19|, |9| and f20fl . 



respectively. For consistency, when including these NLO corrections, we have to use NLO 
unpolarized parton densities and photon fragmentation functions. For the latter we use 
those of GRV [[HJ throughout. The specific choice for the protonic parton densities turns 
out to be rather immaterial, all modern parametrizations leading to essentially the same 
results. The sets of spin-dependent NLO parton densities we want to use have in each case 
been set up in relation to some underlying 'reference' set of NLO unpolarized densities: 
For the distributions of [|J this has been the GRV set [[TJ], whereas the MRS(A') densities 



12 1 have been adopted in ||. For definiteness, we will choose the GRV NLO parton 



distributions |17j as our 'standard' for the unpolarized case, but use the MRS(A') set [12 



when employing the polarized NLO densities of ||. This also implies using the values 
for A^g- (/ being the number of flavors) as implemented in the respective NLO sets, e.g. 
A^r 4) = 200 MeV for the GRV and A^T 4) = 231 MeV for the MRS(A') set, in the NLO 
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expression for the strong coupling ct s : 



*s(Q 2 ) „ 1 A InlnQVAk (4) 
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4vr PolnQVA^ ^o 3 ( ln g2 /A 2_y 
where /3q = 11— 2//3, /?i = 102 — 38//3. As provided for in the sets of flT^ , |L^1 , the number 



of flavors increases when crossing a heavy flavor threshold, and the value for A^i changes 
as a result of the continuity of a s across the threshold. We follow this prescription also for 
the explicit / appearing in the NLO subprocess cross sections. We neglect, however, the 
genuine charm (and, of course, bottom) contributions to the polarized and unpolarized 
cross sections which are tiny at \fs = 39 GeV 113 . 



In the polarized case, next-to- leading order QCD corrections have been calculated for 
the spin-dependent direct subprocess cross sections dAa^ b || [|. We emphasize that care 
has to be taken when combining the NLO sets of polarized parton distributions of 0. |J 
with the NLO expressions for the dAaJ b of f| [| to avoid a mismatch in the factorization 
schemes used. The two-loop splitting functions of 0], which were employed in J|, |(|, 
have been calculated in dimensional regularization in the conventional MS scheme. To be 
more precise, use of dimensional regularization in such a calculation implies to choose a 
prescription for dealing with the Dirac matrix 75 and the Levi-Civita tensor e^ vpcT which 
enter as projectors onto definite helicity states of the involved particles. In |TJ the 'reading 



point' method of |22| with a fully anticommuting 75 was chosen, whereas adopted 
the original definition for 75 of |23|] (HVBM scheme) which is widely considered to be 
the most consistent prescription. It turned out that both calculations |T], [J arrived at 
the same final result for the polarized two-loop splitting functions. The polarized NLO 
parton distributions of M, |6| therefore refer to the conventional MS factorization scheme 



in combination with the HVBM prescription [£3|] (or the one of p2|) for 75, and the 



NLO dAaJb have to be known in the same scheme to make it sensible to use them in 
conjunction with the NLO partons of [B, In fact, the calculation of of the NLO 
corrections to the dAaJ b has been performed using the HVBM prescription, which makes 
the results suitable for our purposes. We note, however, that in a slight deviation 
from the MS scheme was made by factorizing certain finite 'collinear' terms, arising from 
the non-fourdimensional parts of the polarized LO splitting functions as calculated in the 



HVBM scheme, into the NLO spin-dependent parton densities along with the collinear 
singularities ('MSp scheme'). It is straightforward to invert this procedure, i.e., to bring 
back the results of to the conventional MS scheme. Details of this are given in the 
appendix. We note that we refrain from using the results of || for the NLO corrections to 
the dAa^b since it is not obvious whether the 75 scheme used in || provides a consistent 
regularization and can be related in any way to the HVBM prescription. 

Unfortunately, the NLO corrections to the polarized fragmentation subprocess cross 
sections dAa l ab (i = q, q, g) are still unknown as mentioned above. For this reason we stick 
to a pure LO calculation for this contribution, using the corresponding spin-dependent 



LO subprocess cross sections of \2A\ along with the respective LO sets of polarized parton 



distributions of [|, §], and the LO expression for a s (as entailed in (ID by dropping the 



(3\ term) with [17j A^ ^ = 200 MeV. For the photon fragmentation functions we use the 



LO set of |21| in this case. 



Fig. 2a shows the NLO predictions for the spin- dependent cross sections dAa / ' dpTdi] 
as functions of the prompt photon's transverse momentum at c.m.s. rapidity 77 = for 
the four different sets of polarized parton distributions. We also display the unpolarized 
NLO cross section. In Fig. 2b we show the asymmetries All corresponding to Fig. 2a. 
As becomes obvious, All depends strongly on the size and shape of Ag even at large 
Pt- When artificially setting Ag(x, up) = one finds that the asymmetry becomes very 
small and negative for all px and all four sets of polarized parton distributions. Thus the 
differences between the results in Fig. 2b are indeed due to the polarized gluon distribution 
employed. The small negative 'offset' in the asymmetry obtained for Ag(x,[i F ) = turns 
out to be be due to the LO annihilation process qq — > jg and NLO corrections involving 
only incoming quarks and is mainly responsible for the fact that the asymmetries for the 
y Ag = input' and the 'fitted Ag' scenarios of || are much closer to each other at large 
Pt than those for the 'fitted Ag 1 and the ' Ag = g input' scenarios. We also note that 
the full (LO) fragmentation contribution to the polarized cross section is positive for all 
parton sets (apart from set C of || at small to medium pt) and not very sensitive to Ag 
for p T > 7 GeV. 
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Figs. 2b, d show the same quantities as Figs. 2a,c, but now as functions of r\ at pr = 6 
GeV. We note that the rapidity range —1.5 < r\ < 1.5 shown is equivalent to laboratory 
angles of 0.9° < 0j a & < 17°, which roughly corresponds to the range expected to be 
accessible in the HERA-N experiment at this pr ¥M ■ Again, the asymmetries in Fig. 2d 



show strong sensitivity to Ag, even becoming slightly better at large |?7|. 

We have included in the asymmetry plots in Figs. 2b, d the expected statistical errors 
SAll at HERA-N which can be estimated from 



5A LL = 0.17/ yja (pb) . (5) 



This relation has been determined in assuming an integrated luminosity of 240 pb -1 
and beam and target polarizations Pb = 0.6, Pt = 0.8. It includes an overall trigger and 
reconstruction efficiency of 50% but no acceptance correction. The error bars in Figs. 2b, d 
have been obtained using our unpolarized NLO da/dpxdr] in (||), integrated over bins of 
At] = 1, Apx = 1 GeV. The bars have been plotted at the weighted centers of the bins. It 
becomes obvious that the asymmetry should be measurable by HERA-N for p? < 7 GeV 
and for almost all accessible rj, and one should be able to distinguish between different 
scenarios for Ag(x,n^) at 0.1 < x < 0.4. 

The remainder of this paper is devoted to an examination of the reliability of the 
predictions of these very positive findings. We will first discuss the actual size of the 
NLO effects that we have included, and then try to assess the main uncertainties in the 
predictions. 



Figs. 3a, c show the ^-factors 



_ d(A)a NLO 



for the unpolarized and the polarized cross sections for our three sets of spin-dependent 
parton distributions of H. From now on, we do not show the corresponding results for 
set C of H to avoid a proliferation of curves. For the i^-factors the LO cross sections 
in d||) have been calculated consistently, i.e. with LO parton distributions and photon 
fragmentation functions and LO a s . The i^-factors turn out to be quite close to unity in 



the experimentally accessible ranges of pt and r\. Only for the 1 Ag = input' scenario is 
the i^-factor much less than unity around r\ — for all px, indicating rather significant 
NLO corrections in this case. Figs. 3b, d compare the asymmetries All in NLO (as already 
shown in Figs. 2b, d) and LO. By comparison of Figs. 3a,c and 3b, d it can be seen that the 
asymmetries are to some extent less influenced by the NLO corrections than the individual 
polarized and unpolarized cross sections. This is again not true for the 1 Ag = input' 
scenario in Fig. 3b, whose asymmetry changes quite a lot for all pt when going from LO 
to NLO. This feature for a scenario with a small Ag was already observed in [[| |TIJ and is 
due to sizeable negative contributions from direct genuine NLO processes like qq' — > jqq'. 

One major uncertainty in our predictions is expected to come from the fragmentation 
contribution since the parton-to-photon fragmentation functions are experimentally un- 
measured so far, even though very sensible theoretical predictions for them are available 
2T| , p5| (see also In our case, the uncertainty is even larger since, as discussed ear- 



lier, we have to stick to a pure LO calculation for the fragmentation contribution in the 
polarized case, rather than including it on a NLO basis as would be required by consis- 
tency. It seems likely, however, that the individual i^-factor for fully inclusive polarized 
fragmentation is close to unity, just as it turns out to be in the unpolarized caseQ Further 
investigation of this issue is needed in future, requiring a calculation of the NLO QCD 
corrections to the spin-dependent parton-parton scattering cross sections. To assess the 
importance of fragmentation, we show in Figs. 4a, c the ratios 

R = d ( A ) a f™g ^ 

d(A)a dir + d(A)a frag 

in NLO (apart, of course, from the polarized dA<jf rag which is LO). It turns out that the 
fragmentation contribution is generally of 0(20%) at ^/s = 39 GeV and in the ranges 
of pt, f] shown (see also for the unpolarized case), but that it is relatively much 



larger for the 1 Ag = input' scenario around r] = 0, where its importance even rises with 

rising px- This surprising result is, however, in accord with our previous finding that the 

*We note that this is no longer the case if the fragmentation contribution to the isolated prompt 
photon cross section, as measured at very high-energy colliders, is considered. Here the unpolarized 
fragmentation piece turns out to be substantially larger at NLO than when calculated in LO (2(| |27| . 
This feature could set a limitation to the accuracy of similar theoretical predictions for polarized prompt 
photon production at RHIC where the introduction of some experimental isolation criterion will almost 
certainly be necessary. 
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fragmentation part is positive and rather independent of Ag at large pp. Thus in the 
'Ag = input' scenario, where the direct contribution to the cross section is mainly due 
to the small and negative qq — ^ 7$ annihilation process, fragmentation plays an important 
role even at large pp since cancellations between the direct and the fragmentation pieces 
occur. Even though it does not seem likely that the finding of general smallness of the 
asymmetry All for this scenario would become invalid if fragmentation could be included 
at NLO, it means that our predictions for the 'Ag = input' scenario are the least certain 
ones since a strong deviation of the individual .fT-factor for polarized fragmentation from 
unity would affect the predictions for this scenario most. 

On the other hand, we find a very reassuring result in this context when examining 
another major source of uncertainty, namely the dependence of the results on the un- 
physical scales fip and Mp. It turns out that the dependence of all polarized cross 
sections on the fragmentation scale Mp alone is already extremely weak, even though 
fragmentation is only included in LO. This feature, which was also found in [|M| for the 
fully inclusive unpolarized cross section (where NLO fragmentation was used) might be 
evidence for indeed very mild NLO corrections to polarized fragmentation. 

In contrast to this, the dependence of the polarized and unpolarized NLO cross sections 
on the renormalization and initial-state factorization scales is very strong. For instance, 
changing the scales to hr = \ip = Mp = 2p T (we include the change of M F here even 
though it has hardly any influence), the unpolarized NLO cross section decreases by > 50% 
almost uniformly over the whole ranges of pp and 77 (see |TJ]] for a closer examination of 



the notorious scale dependence of the unpolarized inclusive prompt photon cross section). 
Nevertheless, the asymmetries All, which will be the quantities actually measured, remain 
quite untouched by scale changes, as can be seen from Figs. 4b, d. This finding is very 
important since it warrants the genuine sensitivity of All to Ag, implying that despite 
the sizeable scale dependence of the cross sections it still seems a reasonable and safe 
procedure to compare theoretical predictions for the asymmetry with future data and to 
extract Ag from such comparisons. 

In conclusion, we have presented a careful NLO analysis of HERA-N's capability 
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to measure the nucleoli's polarized gluon distribution Ag in inclusive prompt photon 
production with polarized beam and target. The corresponding double spin asymmetry 
An shows strong sensitivity to Ag, and its measurement as well as a distinction between 
various possible scenarios for Ag(x,fi 2 F ) in the range 0.1 < x < 0.4 should be possible 
experimentally. We have also assessed the uncertainties in the theoretical predictions for 
All which appear to be under control. 
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Appendix 



In this appendix we list the changes in the NLO corrections to the direct part of polarized 
prompt photon production which are to be made to transform the results of |J from 
the 'MSp scheme' back to the conventional MS scheme. They are entirely due to terms 
~ e(l — x) in the polarized n = 4 — 2e dimensional LO splitting functions AP^\x) 
(i,j = q,g) as calculated in the HVBM scheme, which were absorbed into the spin- 
dependent parton densities in ||. We emphasize that we must not undo this procedure 
of U for the case of AP$ (x) since in this case the subtraction of the adjacent terms 
~ e(l — x) has also been performed in the calculation of the spin-dependent NLO splitting 
functions where it is demanded by the conservation of the axial non-singlet current [[1], |2|] . 
We therefore only have to reintroduce the effects of the terms ~ e(l — x) in AP^\x) 
7^ {<?<?}) ( as listed in ||) into the polarized NLO cross sections. As a consequence of 
this, only the NLO cross sections for some subprocesses need to be changed. Furthermore 
only the coefficients AC13 of ||, corresponding to the terms involving no logarithms and 
no distributions, are affected. The following terms have to be added to the coefficients 
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AC13 in |J to transform the corresponding polarized NLO subprocess cross sections back 
to the conventional MS scheme: 



qg -» iqg • 



gg -> 7gg : 



qq -> 7gg 



gg' — > 7gg' 



f 2 (l - W) f 2 2 C F , 

— I 2(1 + X)v w + — (1 - 2Xvw t 



-L(2 + Y)v 2 (l-w) 



■(1 - w) ((1 - 2X)(1 + - t; 2 (Xt; 2 + v 3 + Vl )w % ) 



N C X* 

-^-(1 - w) ((1 - 2X)(1 + v)v\ - v 2 (Xv 2 + v 3 + Vl )w 2 ) 



2 Cf / 1 . \ 4/i \ 2 '2^, 



-e^^(l + A> 4 (1 - «,)t^ - e/^-(l + - W ) , 

where = 4/3, X<7 = 3 and the variables v, V\, w, X, F are as defined in |§. For the 
last process, e q and e' q denote the charges of q and q', respectively. The coefficients Aci 3 
for the processes gg — > 'ygg and gg — > jq'q' remain unchanged. 
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Figure Captions 



Fig.l Gluon distributions at Q 2 = 10 GeV 2 of the four NLO sets of polarized parton 
distributions used in this paper. The dotted line refers to set C of 0, whereas the 
other distributions are taken from Q as described in the text. For comparison we 



also show the unpolarized NLO gluon distribution of [1 1 



Fig. 2 a: p^-dependence of the NLO polarized cross section for inclusive prompt photon 
production at HERA-N for the four sets of spin-dependent parton distributions. 
The results are presented at c.m.s. rapidity r\ = 0. The cross section for set C of || 
has been multiplied by —1. For comparison we also show the unpolarized NLO cross 
section. The scales have been chosen to be fiR = fip = Mp = Pt/2. b: Asymmetries 
A LL corresponding to a (here the result for set 'GS C || is shown with its actual 
sign). The expected statistical errors indicated by the bars have been calculated 
according to Eq. (|5|) and as explained in the text. c,d: Same as a,b, but for the 
^-dependence at pp = 6 GeV. 

Fig. 3 a: ^-dependence of the X-factors (as defined in Eq. (§)) for the unpolarized and 
the polarized cross sections at rj — 0. The polarized results are shown for the three 
sets of parton distributions taken from [|J (see text); line drawings are as in Fig. 2. 
b: Comparison of the LO and NLO asymmetries All for the three sets of polarized 
parton distributions taken from |Q (see text). The NLO curves are as already shown 
in Fig. 2b. c,d: Same as a,b, but for the //-dependence at pp = 6 GeV. 

Fig. 4 a: p<r- dependence of the ratios R (as defined in Eq. (^)) for the unpolarized and 
polarized cross sections at rj — 0. The polarized results are shown for the three sets 
of parton distributions taken from || (see text); line drawings are as in Fig. 2. b: 
The scale dependence of the NLO asymmetries All for the three sets of polarized 
parton distributions taken from [|J (see text). The curves for the scales fiR = ftp = 
Mp = pt/2 are as already shown in Fig. 2b and are displayed in the same line 
drawings. For each set the dotted (long-dashed) line corresponds to the scales pp 
(2p T ). c,d: Same as a,b, but for the ^-dependence at pp = 6 GeV. 
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